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ABSTRACT 

Wc present a precise estimate of the bulk virial scaling relation of hales formed via hierarchical 
clustering in an ensemble of simulated cold dark matter cosmologies. The result is insensitive to 
cosmological parameters, the presence of a trace, dissipationless gas component, and numerical res- 
olution down to a limit of ^ 1000 particles. The dark matter velocity dispersion scales with total 
mass as log(crDM(M, z)) = log(1082.9 ± 4.0 km s"!) + (0.3361 ± 0.0026) log(/i(z)M2oo/10^^ Mg), with 
h{z) the dimensionless Hubble parameter. At fixed mass, the velocity dispersion likelihood is nearly 
log-normal, with scatter (Tincr = 0.0426 ± 0.015, except for a tail to higher dispersions containing 10% 
of the population that are merger transients. We combine this relation with the halo mass function 
in ACDM models, and show that a low normalization condition, S'g = erg 3)° '^^ = 0.69, favored 
by recent WMAP and SDSS analysis requires that galaxy and gas specific energies in rich clusters be 
50% larger than that of the underlying dark matter. Such large energetic biases are in conflict with 
the current generation of direct simulations of cluster formation. A higher normalization, 5*8 = 0.80, 

alleviates this tension and implies that the hot gas fraction within rsoo is (0.71 ± 0.09)hjQ^^fli,/Qm, 
a value consistent with recent Sunyaev-Zel'dovich observations. 

Subject headings: cosmology: theory — dark matter — gravitation; clusters: general — intergalactic 
medium — cosmology 



1. INTRODUCTION 

The emergence of the cosmic web of large-scale struc- 
ture is a dynamic, hierarchical process. Galactic-scale 
halos, with formation redshifts greater than one, exhibit 
dynamical maturity through the equilibrium nature of 
the galaxies housed within them. The well-known virial 
scaling relations — TuUy-Fisher for spirals and funda- 
mental plane for ellipticals — are manifestations of ad- 
vanced dynamical age. The largest halos in the universe, 
those that harbor groups and clusters of galaxies, as- 
semble at recent epochs, and for this reason are widely 
thought to be much further from dynamical equilibrium. 

This point of view is supported by manifold obser- 
vational evidence for ongoi ng cluster mergers . Popu- 
lar examples, such as A754 (jHenrv e t al."2004') a nd the 
well-k nown "bullet cluster" 1E0657-56 ( Markevit ch et all 
l2002t l , show multi ple peaks in X-ray emission and galaxy 
number density (jZabludoff fc Zaritskvl Il995f ). In the 
case of the bullet cluster, weak gravitational lensing 
reveals multiple peaks in dark matter that trace well 
the collisionless galaxi es but not the hot intracluster 
gas (jClowe et al.l l2006l ). Similar features are exhibited 
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in beautiful Hubble ACS and Chandra image s of the 
bimodal cluster CL 0152-1357 (|Jee et al.l 120051 ) . where 
a slight lag in X-ray peaks relative to the dark mat- 
ter and galaxies is seen. These signatures, along with 
spatial temperature variations of the sort compiled by 
iBelsole et al.l (|2005[ ). are consistent with those antici- 
pated b y early gas dynamical simulation s of merging 
systems (lEvrard 1990 ; Navar ro et al. 19951 lEvrard et aLl 
[TggllBrvan fc Normaniriggi . 

On the other hand, empirical evidence supports a 
seemingly contradictory point of view of clusters as a 
population of dynamically relaxed systems. Local, X-ray 
flux-limited samples display tight correlations between 
intracluster gas temperature and a number of other ob- 
serva ble features, such as luminosity after c ore exci- 
sion (lMarkevitchlll998t lArnaud fc Evrardl ll999D. isopho- 
tal size (" Mohr fc Evrardl 119971) . intracluster gas mass 
(|Mohr et al.| ll999() . thermal Sunyaev-Zel'dovich decre- 
ment (|Nagail l2006fl. and galactic content meas ured at 
near- infrare d ITn et al.ll2004l: iLin fc Mohij |2004D or red 
bandpasses (jPopesso et al.ll2004D . The level of intrin- 
sic scatter i n these relations i s a few tens of per- 
cent or less (jO'Hara et al.l l2006| ) , making these signa- 
tures, and derived quantities such as gas mass fraction 
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(|Allen et all l2004l : iLaRoque etaOi |2006[ ). a found ation 
for likelihood estimates of cosmological parameters (jVoitl 
[2OOI iLima fcHul[2005l) . 

These lines of evidence paint clusters as a family of 
nearly self-similar, dynamically relaxed systems that can 
be effectively rank-ordered b y a single pa rameter, taken 
by theorists to be total mass (lKaisedll986D . The assump- 
tion of virial equilibrium, which iPeeblesI (|1970[ ) showed 
is rapidly approached after violent collapse and which is 
further validated in this paper, allows gas temperature or 
galaxy velocity dispersion to serve as observable proxies 
for total mass. 

Recent weak lensing studies, however, question the 
regularity of the population by offering evidence for 
large scatter in the lensing mass-temperature rela- 
tion. In particular, the subset of morphologically com- 
plex clusters f ollows a relation o f fset from those of re- 
laxed objects (ISmith et all l2003t ICvpriano et"aI1 l2004l : 
iPedersen fc Dahld |20060 ~ These studies pose the ques- 
tion of whether the massive cluster population should be 
considered a bimodal family. 

The first section of this paper addresses the fundamen- 
tal dichotomy posed above by examining the virial rela- 
tion of massive halos formed in cosmological simulations. 
To the essential question, "Are clusters better described 
as a structurally regular, one-parameter population, or as 
a significantly bimodal combination of 'relaxed' and 'un- 
relaxed' systems?" , the evidence from our simulations 
strongly supports the former, with a caveat that ~ 10% 
of the population are strongly interacting systems. 

All halo finding schemes are challenged by cases in 
which two or more halos are in the process of merg- 
ing into one (| White! 120021 ). Our approach is to use a 
consistent, spherical overdensity definition of dark mat- 
ter halos across a set of computational samples, and di- 
vide the population into primary and satellite halos, with 
the latter being smaller members of spatially overlapping 
pairs. Wc show that the population of halos more mas- 
sive than 10^^ h^^ Mq is dominated by primary systems 
with a regular virial relation. The conditional likelihood 
pio'DulM, z) is very close to log-normal, with only 4% 
dispersion in In ctdm about a robust mean power-law scal- 
ing in mass. Satellite halos are a 5-10% minority of more 
dynamically active objects, with a virial relation offset to 
higher velocities and broader than the primary sample. 

Independent samples from six cosmological codes give 
consistent characterizations of the primary virial rela- 
tion, a testament to the robustness of existing computa- 
tional algorithms. Combining them provides a percent- 
level determination of the slope and intercept, and a 
5% statistical uncertainty in the log-normal dispersion. 
The relation is insensitive to cosmological model, epoch 
and the presence of a minority content of dissipationless 
baryon in halos. 

In the second part of this paper, we demonstrate the 
utility of this calibration by combining it with the mass 
function to predict the space density of dark halos as a 
function of their mean specific energy. Comparing with 
observations of the space density as a function of intr- 
acluster gas temperature and galaxy velocity dispersion 
reveals a tension between expectations from current sim- 
ulations of galaxy and cluster formation and low normal- 
ization cosmologies. 

In §[21 we present our simulation ensemble and method 



of analysis. The ensemble consists mainly of A^-body 
models of dark matter clustering, which provide large 
statistical samples of high mass halos, but we also in- 
clude some simulations that follow the coupled evolution 
of dark matter and coUisional baryons, as well as two 
series of runs that investigate effects of numerical reso- 
lution. Virial scaling relations are analyzed in § [31 in- 
cluding tests of numerical resolution, redshift evolution 
and halo mass definition. Wc explore implications for 
cosmology and cluster physics in § [31 

2. COMPUTATIONAL SAMPLE AND ANALYSIS 

When virial equilibrium is satisfied, the specific ther- 
mal energy of dark matter in a halo of mass M and ra- 
dius R will scale with its potential energy, GM/R. When 
mass is defined using a mean interior density condition, 
then R oc AI^^^ and the kinetic energy scales as M^/'^. 
In keeping with tradition (|Yahil fc VidallfT977l ) . we use a 
one-dimensional velocity dispersion ctdm to express the 
specific thermal energy in dark matter, 

^ 3 

<JUM^ = ^ (1) 
P i=l j=l 

where Vij is the j'^ component of the physical velocity of 
halo member i, the index i ranges over the Np halo mem- 
bers, and V is the mean halo velocity constructed from 
the same Np members. In this work, halo membership 
and total mass are defined using a spherical overdensity 
approach discussed below in § 12.21 

Instead of a classical interpretation of virial equilib- 
rium, expressed as kinetic to potential energy ratio, 
T/W = 0.5, or a boundary-correct ed version thereof 
(|Cole fc Lacevlll996l : [Shaw et al.l2006h . we stress the util- 
ity of the virial theorem as a dimensional tool to con- 
nect total mass to specific energies. The dark matter 
velocity dispersion connects to observable baryonic sig- 
natures, namely the galaxy velocity dispersion (Tgai and 
ICM X-ray temperature Tx, through dimcnsionless ra- 
tios presented in § HI 

2.1. Simulation Ensemble 

We use an ensemble of simulations, consisting of A^- 
body and gas dynamic models of vacuum energy uni- 
verses, listed in Table [Tl along with A^-body realizations 
of other world models given in Table [2l 

The ACDM simulations of Table [1] employ concordance 
parameters, with matter density ilm = 0.3, vacuum en- 
ergy density Q\ = 0.7, and primordial spectral index 
Us = I. Values of the Hubble constant and spectrum 
normalization erg vary across the ensemble, but the virial 
relation presented below is insensitive to these parame- 
ters. References are listed for previously published sim- 
ulations, others arc new to this work. 

The models represent the output of six independent 
simulation codes. Gravitational forces are computed us- 
ing tree algorithms (HOT, PKDGRAV), tree plus mesh 
(GADGET, TREEPM), or particle-particle, particle- 
mesh techniques (HYDRA, P3MSPH). HYDRA, HOT 
and GA DGET are inclu d ed in the code comparison 
study of iHeitmann eFalj (|2005f) . and HOT, HYDRA 
and P3MSPH are among the cod es benchmarked in the 
"Santa Barbara cluster" study of iFrenk et al.l (|1999f ). 
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Sample Np 



L(/i-i Mpc) £(/i-ikpc) Gas? Code 



Remarks /Reference 



HV 


1000^ 


3000 


100 


no 


HYDRA 


2 = and four sky survey outputs 
1 Evrard pt al 1 1120021") 


MS 


2160^ 


500 


5 


no 


GADGET 


Springel et al. (2005) 


MWa 


1024^ 


500 


18 


no 


TREEPM 


White (2002) 


J98 


256^ 


239.5 


36 


no 


HYDRA 


Tenkins et al. (19981 


HOTa 


256^ 


768 


100 


no 


HOT 


new to this work 


HOTb 


256^ 


384 


50 


no 


HOT 


new to this work 


MB 


256^ 


200 


30 


no 


GADGET 


evolved to a = 100 Busha et al. (2006) 


CP 


0.3 - 1 X 10** 


32.5 - 479 


0.06 - 5 


no 


GADGET/PKDGRAV 


23 resimulations Navarro et al. (2004) 


JD 
JB 

MWb 


10^ - 10® 
2 X 1923>'' 
2 X 192^ 


213 

80 - 140 
150 


1.3 - 355 

20-40 

20 


no 

10% 

13% 


PKDGRAV 

P3MSPH 

TREEPM 


resolution series. Dartlv new. Diemand et al. (2004a) 
68 resimulations. Bialek & Evrard (2007) 
White et al. (2002) 


RTM 


2 X 256^'" 


479 


5 


10% 


GADGET 


16 resimulated clusters. Rasia et al. (2004) 



Note. 



Effective particle number in high resolution zone. 



Most of the runs are large-volume representations of 
random regions, but the CP, JD. JB and RTM entries 
are resimulations that focus on a single, dominant halo 
in each run. There are multiple realizations in these sets, 
producing halos that typically span a factor of ten in 
final mass. The exception is the JD resolution series, 
discussed in ij3.2l below. that models the evolution of the 
same initial density field with numerical resolution that 
varies by many orders of magnitude. Overall, the models 
span a factor of one million in particle mass {pL^/Np). 

Although the majority of runs are single fluid and col- 
lisionlcss, the JB, MWb and RTM models follow the evo- 
lution of coupled dark matter and a minority (10 — 13%) 
baryonic component. The physical treatment of the 
baryons in these models is relatively simple, with heating 
due to shocks but no cooling or star formation processes. 

The non-ACDM simulations of Table [2] include two 
simulations of an Einsten-dcSittcr cosmology, Clm ~ 1, 
an open model with = 0.3 and four realizations of a 
dark energy universe with equation of state parameter 
w — —0.8; two each with primordial spectra — 1 and 
0.95. 

2.2. Halo Convention 

We use halo lists identified by individual simulators us- 
ing a common spherical overdensity (SO) halo definition. 
The approach is motivated by a desire to locate hydro- 
static regions centered on peaks in the non-linear density 
field. The individual methods differ slightly in their ap- 
proach to centering and to handling spatial overlap of 
halos, but all employ a consistent definition of mass. 

There are slight differences among the models in the 
treatment of halo centers. In some cases, halo centers 
are identified using a Lagrangian-smoothed estimate of 
the local density, derived from the distance to the 7V* ^ 
nearest neighbor, with 7V~10 (jCasertano fc Hu"3ll985f ). 
In others, a friends-of-friends (FOF) group finder is run, 
and the most-bound particle defines the halo center. 

The mass Ma is defined as the total mass within a 
sphere with mean interior density 3MA/47rrA'^ = Apc{z), 
where Pc{z) the critical density at redshift z. Several 
values of A are in common use, a reflection of the con- 
fusion caused by continual accretio n and mergers ont o 
halos during hierarchical clustering (|Busha et al.ll2005h . 

For most analysis, we employ a fixed critical threshold 



A ==200. In this structurally identical set of halos 

will have virial scaling 

CTDM (X [/l(z)M200]'/', (2) 

with h{z) = iJ(z)/100 km s~^ Mpc^^ the normalized 
Hubble parameter at redshift z. 

We also show results at two other commonly used 
scales with different evolutionary factors. Define Af2oo6 
as the mass within a sphere encompassing a mean mass 
density of 200 times the background matter density 
Pm{z). For this background mass, one expects a scal- 
ing of the form 

(TDM cx [(l + z)3/2M2oofc]'/'. (3) 
Finally , analytic solutions for spherical perturbation evo- 
lution (jLahav et al.lll99l[ l motivate the use of a variable 
threshold Av ar(2:) against the critical density. We use 



Eke et all (|1996f) . but scale to an asymp- 



the values of 

totic value of 200, rather than 178 at high redshift (the 
Qjn ~ 1 limit). This implies a present-epoch value 
Avar(O) = 115 in a concordance cosmology. Defining 
fvar{z) = Avar (2)/ Avar (0), the expcctcd rcdshift scaling 
for the virial relation is 

11/3 



"■dm 



cx [fll,'{z)h{z)M,.. 



(4) 



Although this scale has typically been referred to as 
the "virial mass" Mvw, we employ the label var to 
avoid misinterpretation of this mass as th e unique scale 
at wh ich virial equilibrium is satisfied. ICole fc LacevI 
(|1996f ) demonstrate that the virial ratio 2T/W is a weak 
function of radius ne ar the scales discussed here, and 
ICohn fc White! (|2005f ) show that it is also a weak func tion 
of redshift in a ACDM cosmology. IShaw et all (|2006( l ex- 
amine the pressure-corrected virial ratio 



/?vir = 



2T-Es 
W 



1, 



(5) 



with T the kinetic energy. W the gravit ational potential 
energ y and Eg a boundary pressure term (|Chandrasekhail 
Il961h . and find that /Jvir is nearly constant in the range 
A'^ 100 — 500. The results presented in § 13. 3 1 confirm that 
tight virial scalings exist across the full range of scales 
discussed above. 

We use corresponding notation for the radial sizes of 
the above mass definitions: r2oo is the radius that defines 
-M200, ?'2006 defines Af2oob, and 7^ 
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Table 2. non-ACDM Af-BODY Simulations 



Sample 




L(/i-i Mpc) 


£(h ^ kpc) 


Code 


Romarks/Rcfcrcnco 


HVt 

J98r 
MWc 

MWd 


1000^ 

256^ 
256^ 

256^ 


2000 

239.5 
256 

200 


100 

36 
24 

30 


HYDRA 

HYDRA 
TREEPM 

TREEPM 


f2m = 1 , z = and four sky survey 
outDuts lEvrard et al.l 1^6%) 
n„ = l. Jenkins et al. fl998) 
^^ = 0.3 

4 realizations, Qm = 0.357, 
w = ~0.8, ns = l,0.95 



2.3. Interacting Halos 

Unlike percolation algorithms, which have the virtue of 
creating exclusive group assignments for simulation par- 
ticles, the SO algorithm potentially allows a particle to 
belong to more than one group. The halo finding schemes 
we employ differ in their treatment of overlapping halos; 
the scheme used with the HV, J98 and HOT simulations 
allows it, while the others do not. In all cases, however, 
an exclusivity condition applies to halo centers; the cen- 
ter of one halo cannot lie inside the spherical boundary 
of another. 

In the cases where overlap is allowed, we label the 
larger member of an overlapping pair the primary halo, 
and the smaller member the satellite. Note that satellites 
should not be confused with sub-halos; the latter are usu- 
ally sub-units lying within a larger halo. The satellites 
here are merely smaller members of pairs whose spherical 
volumes intersect. Figure 1 shows the fraction of satel- 
lite halos identified in the ACDM HV model, the simu- 
lation with the largest population of massive halos. At 
z = 0, satellites account for 12% of the population above 
10^* Mq. The fraction dechncs roughly as M~^/^ at 
higher masses. Although more massive systems are more 
strongly clustered, they are also much rarer, and the lat- 
ter dominates to make the satellite frequency a decreas- 
ing function of mass. Somewhat surprisingly, three out 
of every hundred Coma-sized {10^^ ^q) halos will 
lie within r2oo of a more massive neighbor. 




Fig. 1. — The fraction of halos classified as satellites 
— smaller members of distinct overlapping halo pairs — 
as a function of mass in the ACDM HV model at z = 
(dashed) and in sky surveys extending to 2; = 0.54 (light) 
and 1.25 (bold). 



Also shown in the figure are satellite fractions for sky 
survey outputs of the HV simulation that cover Stt stera- 
dian of sky to z = 0.54 and tt steradian to z = 1.25 
(|Evrard et al.|[200^ . At higher redshift, the more dilute 
high mass halos are less likely to overlap. The satellite 
fraction above 10^^ Mq is 8% for the sky survey ex- 
tending to z = 1.25. 

Our estimate of 10% for the fraction of strongly in- 
teracting haJos_Js_sornewlm higher than the 3.4% de- 
rived by I Shaw et al.l (|2006[ ) in their analysis of a com- 
putational sample of 2000 well-resolved (> 10^ particle) 
ACDM halos more massive than 3 x 10^^ h'^ Mq. From 
their ma in halo sample, ide ntified with a percolation al- 
gorithm (jWeller et al.ll2005f l tuned to A~50, they define 
a dynamically disturbed, or irregular, subset using a con- 
dition on the virial ratio, /?vir < —0.2. Shaw et al. note 
that the distribution of /3vir is continuous and that their 
choice of threshold is motivated by a desire to flag ex- 
treme cases. From their Fig. 1, it appears that raising 
the /?vir threshold to a value —0.15 would roughly double 
the fraction of irregular halos, bringing it more into line 
with our estimates based on spatial overlap. 

3. THE VIRIAL RELATION OF DM HALOS 

We show that massive dark matter halos adhere to a 
virial scaling relation of the form 

aDM(M,z) = crDM,15 I -^^15 I , (6) 

with (TDM. 15 the normalization at mass 10^^ Mq and 
a the logarithmic slope. 

The final result of this section, derived by combining 
all well-resolved, ACDM samples of primary halos, is 
a percent-level precision estimate of the fit parameters, 
CTDM,i5 = 1082.9 ± 4.0 km s'^ and a = 0.3361 ± 0.0026. 
Tabic [3] lists these parameters along with the standard 
deviation of In ctdm about the best fit relation. We dis- 
cuss this result in § 13.51 

We begin by establishing agreement among codes for 
the present-epoch relation in both ACDM and non- 
ACDM cosmologies, then move on to test numerical con- 
vergence using the JD and CP samples. Time evolution 
and scale dependence of the virial relation are presented 
in § 13.31 followed by a brief examination of its future be- 
havior in the deSitter phase of a ACDM cosmology, the 
era of genuine halo equilibrium. 

3.1. Present-epoch Relation 

Figure [5] shows the virial scalings of the halos in the 
ensemble of ACDM runs at z = 0. Each panel shows 
the internal velocity dispersion, equation (1), for primary 
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Fig. 2. — The dark matter halo virial relation for ACDM models, defined in Table[T] at 2 = 0. Point styles are primary halos (dots) and 
satellites (open circles). One thousand, randomly chosen halos are shown for samples larger than this size. The top two rows are A^-body 
simulations of dark matter, the bottom row includes gas as a second, coUisional fluid. The dashed line in each panel is the least-squares fit 
for that model (Table [6ll. The solid line in all panels is the global best-fit with parameters listed in Table [3] 



Table 3. A = 200 Primary Halo Virial Relation (A^is > lO") 



Parameter Value 

slope, a 0.3361 ± 0.0026 

intercept, ctjjm.is 1082.9 ±4.0 km s^^ 

scatter, (gf,,^)^'/^ 0.0426 ili 0.0015 



halos with mass A/200 > \-^^^h-^ M©. In addition, the 
HOT, J98 and HV panels also display values for satellite 
halos, shown as open circles. For clarity, samples with 
more than 1000 halos are subsampled to that level. 

For each model, we perform a least-squares fit to 
log((TDM) = log(crDM,15) + a log(/i(z) Af2oo/10^^ Mq) to 
all primary halos above the 10^"^ h'^ Mq limit. Although 
using equal weight per halo puts more emphasis on sys- 
tems near the mass cutoff, we find that fits to binned 
versions of the data produce similar results. The least 
squares fit for each model is shown by the dashed line 
in each panel, while the solid line shows the global fit of 
Table El 

Best-fit parameters for the primary halo samples are 
given in Table [6l Also listed is A'haio, the number of 
halos above the 10^'' Mq limiting mass, and mass 
resolution A^is, the number of particles in a 10^^ /i"^ Mq 
halo. The listed uncertainties are \a statistical errors 
from the least-squares fits. 

The ensemble of models covers a wide range in sam- 
ple size and mass resolution. The Hubble Volume sam- 



ples contain one-half million halos, but the minimum 
mass limit is resolved by just 45 particles. At the 
other extreme i s the M illennium Simulation sample (MS, 
ISpringel et all (j2005l) ) of nearly 2000 halos with 2500 
times better mass resolution than HV. In all cases, the 
primary halo population adheres to a tight relation, with 
rms deviation of 4 — 5% in (Tdm at fixed mass. Indepen- 
dent codes produce consistent results for this form. 

Models with poor resolution exhibit significantly 
steeper slopes than those at higher resolution. For exam- 
ple, a = 0.357 ± 0.001 in the HOTa run with N15 = 450, 
while a ^ 0.335 ± 0.002 in the MWa model with 200 
times better mass resolution. We show in § 13.21 that this 
steepening is consistent with the systematic error intro- 
duced by numerical resolution at the low mass limits of 
the HOTa and HV models. On the other hand, our reso- 
lution tests below support the finding that these models 
have sufficient resolution to accurately measure the in- 
tercept at 10^^ h-^ Mq. 

Figure [3] plots the conditional likelihood, p{aDM\M), 
determined from the residuals about the best-fit power- 
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Fig. 3. — Distributions of deviations in IncrjjM about the best-fit relations shown in Fig. [2] for primary (open histograms) and satellite 
halos (shaded histograms). The solid curve in each panel is a normal distribution with scatter derived from the primary halos of each 
model and given in Table [6] The normalization of the satellite distribution is enlarged by a factor of 5 for clarity. A minimum mass of 
10^^ Mq is applied, except for the low resolution models HV and HOTa, for which 5 X lO'^* Mq is used. 



law relations for the samples of Figure [2l For the low 
resolution runs HV and HOTa, we increase the mass limit 
in this figure to 5 x 10^^ Mq; all other panels remain 
limited at lO^/i"^ Mq. 

The velocity dispersion likelihood for primary halos is 
nearly log-normal with standard deviation (5f„o.)^^^ ^ 
0.04. There are clear signs of non-zero skewness, par- 
ticularly in the MS and MWa samples that have the 
largest number of well-resolved halos. The median val- 
ues of In(crDM) are 0.08 and 0.10 times the standard de- 
viations of these samples, and the skewness, N'^(x — 
- ^)^)^/^ is 0.74 ± 0. 12 and 0.85 ± 08, re - 
spectively. The analytic model of lAfshordi fc CenI ()2002f ) 
anticipates skewness, but it also predicts that the vari- 
ance should decline with increasing mass, a feature that 
is not seen in our data. 

The likelihood p(crDM|Af) for satellite halos in the HV, 
HOT and J98 models is displaced to higher values, by 
0.13 in mean ln((TDM), and is broader than the primary 
population by a factor two. These features are consistent 
with a picture of satellites as merging systems. Since 
the simple density threshold definition used for halos 
does not take particle binding energy into account, the 
satellite systems typically consist of a mixture of some 
locally bound material along with a hotter component 
contributed by the neighboring primary. More detailed 
group-finding treatments, such as those that combine a 
percolation algorit hm with dynamical identificati on of 
sub-halos (Sprinac l et al.ll2001l : IWeller et al.ll2005f ). will 
be required to clarify the nature of satellites. 



The virial relations of non-ACDM models are shown 
in Figure 4. As in Figure [21 dashed lines are best fits to 
each model and solid lines give the global ACDM relation 
of Table [3l Despite the fact that the assembly history of 
halos is sensitive to cosmology, the resultant virial rela- 
tions are remarkably robust. The power-law fits of these 
models are consistent with that of the ACDM family to 
within a few percent in intercept and slope. 
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Fig. 4. — The z = virial relations of the non-ACDM 
models listed in Table [2l Point and line styles are the 
same as Figure [2l Fit parameters for primary halo sam- 
ples are given in Table [G] 
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3.2. Numerical Resolution 

Simulations of random cosmic volumes have the benefit 
of producing statistical samples of halos that are unbi- 
ased in the large volume limit, and can be made so by 
appropriate mode-sampling of finite-volume realizations 
(|Sirkall20"05l ) . The principal drawback of this approach is 
the differential mass resolution imposed by the fixed par- 
ticle mass; smaller halos are simply composed of fewer 
particles. The HV and HOTa samples push numerical 
resolution toward its lowest practical limit {Np > 45). 

On the other hand, constrained initial condi- 
tions ([Bertschingeij 119871) and resimulation techniques 

avarro fc Whitd Il994l : iTormen eTall 119971) can pro- 
duce ensembles with fixed resolution across a range of 
mass scales, but halo sampling under this approach need 
not be uniform or volume complete. 

We find that the resimulation and large volume sam- 
ples produce consistent normalizations for the virial re- 
lation. However, estimates of the slope show sensitivity 
to numerical resolution. The poorly-resolved models, HV 
and HOTa, have a ~ 0.36, equivalent to a ~7% reduction 
in velocity dispersion at 10^'' Mq with fixed ctdm.is 
normalization. 

To address whether this reduction in velocity disper- 
sion is to be expected purely from numerical resolution, 
we turn to the JD ensemble, a resimulation series o f a 
single halo originally used bv lDiemand et al] (|2004a[ ) in 
a study of inner density profiles. The published work 
contained high resolution resimulations. We include here 
lower resolution realizations that use the same initial per- 
turbations truncated at progressively lower Nyquist fre- 
quencies. 

Table 2] lists the halo properties of this series. The 
12 runs span nearly six decades in mass resolution, with 
the number of particles within r2oo, ranging from 
33 to 10''. At lower resolution, the loss of total power in 
the realized density field delays halo collapse, leading to 
lower values of the mass M200 at z = 0. The velocity dis- 
persion also tends to decrease as the resolution degrades. 
The exception is the realization with iV2oo = 82. From 
visual inspection of the output, this system differs from 
the others in terms of its dynamical phase. It is experi- 
encing a major merger while systems in the remainder of 
the series are relatively quiescent. 

For each member of the series, we employ M^/^ scaling 
to determine the effective value of the intercept ctdm.is- 
The results are listed in Table |4] and are plotted against 
numerical resolution 7V200 in Figure 5. The inferred ve- 
locity dispersion normalization is extremely robust. It 
lies within ±2% of the highest-resolution value until the 
number of particles within r2oo falls below ^ 500 parti- 
cles. At the level of 30 particles within r2oo, the drop 
in ctdm.is has grown to 15%, still a modest reduction 
given the extremely degraded resolution. 

For the HV and HOTa simulations, these results sug- 
gest an underestimate, at the ~ 10% level, in velocity 
dispersion at 10^* Mq, or 45 particles. At the nor- 
malization mass scale of IQ^^ Mq, the inferred bias 
should be closer to ~ 1%. The differential effect of this 
bias is consistent with the degree needed to tilt the HV 
and HOTa relations by ~0.03 in slope, while leaving the 
intercept ctdm.is little affected. 

Determination of the slope a is best done with high- 



resolution simulations that cover a wide dynamic range 
in mass. The CP series satisfies this design criterion, 
covering five orders of magnitude in mass with million- 
particle per halo simulations. 

The virial relation for the CP series, shown in Fig- 
ure 6, has best-fit slope a — 0.331 ± 0.002, a value con- 
sistent with the self-similar expectation of 1/3. Given 
the well-documented evidence for departures from self- 
similarity across the mass spectrum, the 1% agreement 
in slope is certainly surprising. Compared to cluster- 
scale halos, those of dwarf galaxies are rounder, more 
central ly concentrated, a nd have more isotropic particle 
orbits (jPrada et al.ll2006l ). Despite these differences, the 
self-similar form of the virial relation is respected across 
the range 10^" to lO^^/i"^ Mq. The intercept and scat- 
ter in ctdm are consistent with the ensemble values (see 
Table [SD. 
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Fig. 5. — Inferred virial relation normalization as a func- 
tion of number of particles within r2oo for the JD resolu- 
tion series of a single halo. The shaded region denotes a 
2% range centered on the highest resolution value. 
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Fig. 6. — The virial scaling for the CP set of 23 very 
high resolution halos spanning dwarf galaxy to cluster 
mass scales. The solid line shows the best-fit power law, 
and the inset shows residuals about the fit. 
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Table 4. JD resolution series 



-'V200 


[ iU fl 




D iVi , i 


AA^ \ ('km s"-'- 




10519395 


0.159 


601 


1108 


4756791 


0.171 


619 


1114 


1386543 


0.168 


608 


1101 


1362015 


0.165 


606 


1104 


158962 


0.154 


581 


1083 


5806 


0.152 


575 


1077 


698 


0.146 


582 


1105 


220 


0.155 


550 


1023 


193 


0.137 


513 


995 


82 


0.137 


573 


1113 


34 


0.112 


461 


956 


33 


0.107 


453 


953 



3.3. Evolution and Mass Measure 

The preceding analysis establishes a well-defined virial 
relation for massive halos at the present-epoch. We now 
address how this relation evolves in time. Along with 
the default M200 measure, we also examine the evolution 
based on the other mass scales discussed in § 12.21 

Figure [7] shows the virial relation for halos sampled in 
the redshift range z = to 1.5. The data include sky 
survey samples from HV models of ACDM and rCDM 
cosmologies, along with samples at 12 discrete redshifts 
from gas dynamic simulations of ACDM clusters. Fits to 
these data are given in Tabic [6l 

In all three cases, the virial relation scales according 
to the expectation of equation ([2]). By virtue of offering 
a low-dispersion estimator of the potential well depth, 
cdm, one might argue that the measure h{z)M2m de- 
serves consideration as the "virial mass" of a halo. 

On the other hand, the regularity of halo structure im- 
plies that this scaling relation is not unique. In Figure [8l 
we compare virial relations for the three mass measures 
discussed in § 12.21 For this exercise, we use a single halo 
sample taken from the one of the deep octant sky survey 
outputs of the ACDM HV simulation (the NO sample) . 
We employ a mass limit 3 x 10^^ M© as a compro- 
mise between sample size and numerical resolution. Fits 
to the data are listed in Table [5l 

The mass measures sample regions of differing radial 
extent, with M200 the innermost and A/2006 the outer- 
most scales. The virial relation intercept drops at larger 
radii, a reflection of our choice of normalization at a fixed 
mass of 10^^ Mq and of the declining velocity disper- 
sion profiles in indivi dual halos beyond the scale radius 
([Navarro et al.l[T"997[ ). The slope a, although biased high 
due to resolution effects discussed above, is consistent for 
all three mass measures. 

The fractional scatter in the likelihood p{(J'dm\M^z) 
grows slowly as the radial scale is increased, from 0.0487 
at M200 to 0.0559 at M2oob- Jackknifc errors in the scat- 
ter are 0.0005, implying that the small change in scatter 
is significant. The longer dynamical times in the outer 
portions of halos is the likely cause of the increasing vari- 
ance. 

The trend of decreasing variance at smaller scales is 
unlikely to continue to the core regions, r < 0.1r2oo- 
At such small radii, the influence of minor mergers and 



sub-halo transits can be large, and the variance is there- 
fore likely to increase. The HV-zO sample at A = 500 
has scatter 0.042, slightly smaller than the 0.0439 value 
at A = 200. This evidence suggests that the minimum 
variance scale the virial relation lies within the A = 500 
critical surface. 



3.4. Asymptotic Behavior in the ACDM case 

As a final comm ent on evolution, we use the MB run 
(jBusha et al.ll2006l ) to explore the behavior of the virial 
relation in the future of a ACDM cosmology. We exam- 
ine the behavior at scale factor a/ao = 100, 53 Gyr 
into the future, and well into the deSitter phase of 
the vacuum-dominated era when linear growth of large- 
scale structure has shut down. The physical scale of 
the cosmic web is exponentially stretched, and the em- 
bedded, non-linear halos evolve in increasing isolation 
(jNagamine fc Loebl |2003[ ). In short, the web has con- 
densed into a spray of isolated halo droplets. 
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Fig. 7. — The halo virial relation at non-zero redshift for the models indicated. Point and line styles are the same as Figure [2] The HV 
models are continuously sampled along light-cones extending to z~1.5 while the JB runs are sampled at twelve epochs, spaced equally in 
time, with z < I. 
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Fig. 8. — Virial relation of primary halos in the HV-NO sky survey sample, with masses defi ned b y (left to right) fixed critical, variable 
and fixed background density thresholds. The appropriate prefactors of mass arc described in ^ 12.21 
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Fig. 9. — The future virial relation of an ACDM cosmol- 
ogy. Points show halos at a~ 100, and the solid line gives 
the best fit power-law relation. The dashed line shows 
the relation at a < 1. The inset gives deviations about 
the best fit (shaded) and the present-epoch (unshaded) 
relations. Loss of unbound material leads to lower veloc- 
ity dispersions and dynamical relaxation decreases the 
scatter in velocity dispersion to below 2%. 



At the present epoch, the ensemble-averaged, radial 
phase-space structure of halos has two zero velocity 
surfaces: an inner hydrostatic radius and an outer 
turnaround radius. Within the wide mixing zone be- 
tween them lie the mass scales discussed in § 12.21 In the 
near future of a ACDM universe, the hydrostatic surface 
grows and merges with the turnaroun d radius, creating a 
single, well-defined cluster edge (Bu sha et al.|[20051 near 
the classical turn-a round radius defined by critical de n- 
sity contrast A~6 l|Gunn fc Gottlll972t lPeebleslll98GD . 

Figure 9 shows that the a = 100 relation is a slightly 
lowered and substantially tighter version of its present- 
day counterpart. The intercept, ctdm.is — (1053 ± 
6) km s^i, is 3% lower than the present-epoch relation 
for this simulation, shown in Figure [2l The lowering re- 
sults from the loss of a small fr action of unbound mass 
that occurs as accretion ceases (jBusha et al.ll2005l ). 

The scatter in p{(Tj:,m\M) is 1.9±0.1 percent at a=100, 
meaning the variance has fallen by more than a factor of 
four compared to the present epoch. The variance today 
must therefore be dominated by transient phenomena as- 
sociated with accretion and mergers, not by long-lived 
structural differences such as variations in halo shape or 
concentration. Such structural differences persist into 
the near future. 
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Table 5. Virial Relation for Alternate Mass Measures 



Mass 




a 


""DM, 15 
(km s-1) 


\ lllCT' 


M200 


4175 


0.352 ±0.003 


1080 ± 2 


0.0487 ± 0.0005 


Mvar 


5308 


0.355 ±0.002 


982 ± 1 


0.0527 ± 0.0005 


-'^2006 


6363 


0.355 ±0.002 


880 ± 1 


0.0559 ± 0.0005 



Note. — Derived from HV-NO sky survey samples, mass-limited at 3 X 10^* h ^ Mq in each measure. 



3.5. Calibration Summary 

The virial scaling parameters for primary halos, shown 
in Figm'es 2, 4, 6 and 7, are summarized in Table [Sj 
Figure 10 shows the dependence of the fit parameters on 
the degree of numerical resolution, characterized by the 
7Vi5, the number of particles within a 10^'' h"^ Mq halo. 
The slope and, to a lesser extent, the scatter are biased 
high in the lowest resolution simulations, but converge for 
A^i5 > lO'^, or 1000 particles at the 10^'^ h'^ Mq limiting 
mass. The intercept ctdm,15 shows little sensitivity to 
resolution. 

The high degree of consistency among the models moti- 
vates us to combine them to produce an overall fit. We do 
so using a bootstrap approach on the ACDM models of 
Table[6]with resolution iVi5 > 10*. We make one adjust- 
ment to the normalization error, based on the recogni- 
tion that computational sources — force anisotropics and 
time integration error in the A'^-body algorithms — will 
produce a fioor in (Tdm,15 uncertainty. For modern algo- 
rithm s, typical errors in kine tic energy are at the level of 
0.5% (jEfstathiou et al.lll985[ ). We therefore add a fixed 
0.25% fractional error in quadrature to the statistical un- 
certainties in (TDM. 15- Our bootstrap approach generates 
a large number of trials, each realizing a random selec- 
tion of fit parameters assuming Gaussian errors. We then 
quote the mean and standard error determined from the 
trial distributions. 

The resultant global values, quoted above in Table [3l 
are shown by the bold lines in Figure 10, and shaded 
regions show 90 percent confidence errors. The slope 
and intercept are determined to better than one percent 
precision; the uncertainty in the intrinsic variance (Sf^^) 
is slightly less than ten percent. 

The global fit values are robust to the choice of Ni^ 
threshold. Including the less well-resolved models in the 
fit raises the slope to 0.3393, while the intercept and 
dispersion change by less than la. To the extents probed 
by the simulation ensemble, the fit values are insensitive 
to cosmological model, epoch, the presence of a trace 
baryonic component, and the computational algorithm 
used in the simulation. 

The uncertainties in the virial fits are driven by sta- 
tistical errors. A concerted effort among computing con- 
sortia to federate large ensembles of moderate- to high- 
resolution simulations could readily push statistical er- 
rors in (TDM, 15 below the level of current algorithmic 
accuracy. This eventuality would drive a need for in- 
proved understanding of energy conservation in existing 
algorithms, potentially stimulating the development of 
higher-accuracy cosmological codes. 

In the next section, we use the percent-level accuracy in 
the virial relation as an anchor for the observable specific 



energies in rich clusters. 

4. IMPLICATIONS 

The virial relation establishes a strong connection be- 
tween mass and specific kinetic energy in dark matter. 
Since experimental detection of dark matter rem ains elu- 
sive even m our own galaxy (jEUis et al.ll2005f ). direct 
measurement of these velocities is currently out of the 
question. However, the dark matter potential well depth 
sets the scale for baryonic physics, and we consider here 
simple paramcterizations for the the X-ray temperature, 
the galaxy velocity dispersion, and the ICM mass frac- 
tion. 

4.1. ICM and Galaxy Components 

The mix of baryonic phases in halos is a complex 
function of time and halo mass, driven by compet- 
ing, non-linear cooling and heating processes. Rich 
clusters have an advantage in that they represent re- 
gions of accelerated ga lactic development (jBowedllQQll : 
iKauffmann et al.l I1993D . Compared to a typical re- 
gion in the field, star and massive black hole forma- 
tion pla y out at earlier cos mic times in proto-clustcr 
regions (jSpringel et al.l 120051 ). This, along with the re- 
duction in cooling rates caused by a growing virial tem- 
perature, allows heating to dominate coolin g early on, re- 
sultin g in a low galaxy formation efficiency (jPearce et al.l 
|1999() . The baryonic content of rich clusters is thus ob- 
serve d to be only 10 - 15% stars and cold (T < 10^ K) 
gas (|White et al.l Il993l : lEtFori ,200a). Recent studies 
that include a diffuse light contribution emphasize the 
mass-dependence of this fraction; stars and hot gas are 
( comparable mass com ponents at M200 = 10^'' Mq 
(jGonzalez et al.ll2007t ). 

The mean halo ICM mass fraction within ta will differ 
from the global baryon fraction rif,/J7m, due partly to 
the mass in cold gas and stars but also to differences in 
the radial structure of the dark matter, ICM and galaxy 
components. We parameterize these effects separately. 

We first introduce the mean halo baryon bias 



(7) 



where Mh is the baryon mass within ta- In addition, we 
introduce the baryonic ICM mass fraction within ta 

MjCM 



/] 



ICM 



(8) 



With these definitions, we can relate the local ICM 
mass fraction to the global baryon fraction 



JiCM = -YT — — ./ICM 1 7^ 

AiA V i 'm 



(9) 
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Fig. 10. — Resolution dependence of the primary halo fit parameters for the ACDM 2 = data without/with gas (filled/open circles), 
ACDM 2 > data (open squares), and non-ACDM models (asterisks). The bold horizontal line and shaded region in each panel are global 
means and 90% confidence uncertainties for the ACDM models determined from bootstrap resampling. 



Table 6. A = 200 Primary Halo Virial Relation Fits 



Sample 


-'Vhalo 


Afl5 


0"DM,15 


a 


\ Incr / 


HV-zO 


23636 


450 


1073.5 ±0.5 


0.354 ±0.001 





0439 ± 0.0002 


MS 


1938 


1162800 


1093 ± 4 


0.340 ± 0.002 





0386 ± 0.0006 


MWa 


4524 


103180 


1086 ± 3 


0.336 ±0.001 





0409 ± 0.0006 


J98 


319 


14660 


1080 ±11 


0.338 ± 0.006 





049 ± 0.003 


HOTa 


427 


450 


1079 ± 3 


0.349 ± 0.007 





045 ± 0.002 


HOTb 


1358 


3400 


1090 ± 5 


0.343 ± 0.002 





045 ± 0.002 


MB 


133 


25190 


1072 ± 13 


0.336 ± 0.007 





036 ± 0.004 


JB 


67 


-15000 


1053 ± 8 


0.316 ±0.007 





037 ± 0.003 


MWb 


158 


-25000 


1075 ± 14 


0.337 ±0.007 





040 ± 0.003 


RTM 


16 


-300000 


1121 ±28 


0.361 ±0.032 





058 ± 0.013 


HVr-zO 


25743 


450 


1079.5 ±0.5 


0.355 ±0.001 





0500 ± 0.0003 


J98t 


652 


4400 


1085 ± 9 


0.337 ± 0.004 





050 ± 0.002 


MWc 


586 


18020 


1102 ±7 


0.341 ± 0.004 





046 ± 0.002 


MWd 


984 


21170 


1101 ±6 


0.342 ± 0.003 





0419 ± 0.0008 


HV-ss 


20553 


450 


1077.5 ±0.7 


0.350 ±0.001 





0472 ± 0.0002 


HVt-ss 


13049 


450 


1079 ± 1 


0.361 ±0.002 





0503 ± 0.0002 


JB-zl 


777 


-15000 


1075 ± 3 


0.330 ± 0.003 





045 ± 0.002 


CP 


23 


-10'^ 


1087 ± 16 


0.331 ± 0.002 





044 ± 0.006 



Note. — Least-squares fits use equal weight per halo above limiting mass 10^* h ^ M0 (5 X 10^* h ^ M© for HOTa and all HV runs). 
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To describe baryon energetics, we introduce dimension- 
less specific energy ratios that enable connections to op- 
tical and X-ray observations. Let 



bv = 



CTDM 



(10) 



be the ratio of galaxy velocity dispersion to that of 
the dark matter, al s o known as the v elocity bias 
(ICarlberg et all 119901 : fevrard et all Il994l: IColfn et al l 



20001: ISorineel et all [20011 iFaltenbacher et all l2005l: 



Faltenbacher fc Diemandll2006 ) . The unbiased case, 



1, is expected if both dark matter and galaxies behave as 
collisionless fluids, and if the galaxies fairly sample the 
mass- weighted phase structure of the dark matter. 

The velocity structure of clusters is sensitive to galaxy 
type, with late-t ypes having roughly 10% higher a^sd 
than early- types (jSodre et al.lll989l : lde Theije fc KatgertI 
Il999f ) and the br ightest galaxies hav ing substantially 
lower dispersions (jBiviano et al.lll992[ ). We do not ad- 
dress these detailed dependencies in this paper, but con- 
sider only the velocity dispersion of magnitude-limited 
samples. 

For the hot gas, we use the specific, thermal energy 
defined by the X-ray temperature Tx to form the ratio 

2 _ kTx/^J.mp 



11 



with /i the mean molecular weight and the proton 
mass. N ote that 6^ is the inverse of the c lassical /3 pa- 
rameter (jCavaliere fc Fusco-Femianolll976t ). That work 
envisioned a potential well dominated by the mass in 
galaxies, so their original measure is now better known 
as the spectroscopic beta parameter 

g = i£l - (12) 



kTx/fJ-mp 



by 



This quantity has the benefit of being directly observable. 
From a heterogeneous sample of 109 clusters, [Xuc & Wul 
((2000h find a mean value /3spcc = 1.0 at kTx = 6 keV. 

4.2. Direct Simulation Constraints 

Direct simulations of the ICM and galaxy components 
provide independent estimates of the above m easures. 
Two signatures exhibited by early simulations (jEvrardl 
119901 : iNavarro et al.l [T995[ ) have remained robust, even 
as resolution has increased and new modeling algorithms 
have emerged. One is that the hot ICM is not completely 
thermalized. Within r'2oo, kinetic pressure in ordered and 
turbulent flows contributes 10 — 15% relative to thermal 
pressure. The second is that, under purely gravitational 
evolution, the local baryon fraction within this scale is 
slightly depleted relative to the global value. 

The 'Santa B arbara Cluster' comparison study 
(|Frenk et al.l |1999[ ) applies twelve gas dynamics codes 
to the study of a single cluster realization. The codes 
agree that the ICM is not completely thermalized in the 
evolving potential. The resultant mean specific thermal 
energy of the gas within r2oo is 



(0.87 ±0.04) {Tx/T^}. 



(13) 



The quoted error is the standard deviation among the 
different codes^, and the term in brackets arises because 

^ Excluding the lowest resolution model t hat had only a few 
hundred volumetric cells within the cluster, see IPrenk et aUllT999l '). 



the simulations employ a mass- weighted temperature Tm 
while the definition of b^ is with respect to spectroscopic 
temperature measured by X-ray observations. We em- 
ploy below a value (Tx/Tm) = 1.10 ± 0.05 measured 
within r5oo by iNagai et al.l (|2006[ ) using detailed mock 
Chandra observations of simulated clusters. 

The Santa Barbara solutions also provid e an estimate 
of the halo baryon bias, T = 0.86 ± 0.06. iPearce et al.l 
(|1994f l offer evidence that the depletion arises from en- 
ergy transfer into the ICM from the dark matter during 
major mergers. Recen t ART and GADGET simulations 
(jKravtsov et al.ll2005l ) find values consistent with this es- 
timate when no gas cooling is allowed, but models with 
cooling and star formation tend to be overabundant in 
baryons, T ~ 1.05 at r^QQ. Note that these simulations 
do not include the effects of AGN feedback. 

The fraction of baryons that condense in to stars and 
cold g as is difficult to model accurately. lEttori et al.l 
(I2006h analyz e the 480^-particle GADGET models of 
IBot gani et al\ (|2004f l plus additional resimulations and 
find that, while the overall baryon fraction within r2oo 
is relatively robust, T = 0.94 ± 0.03, the hot gas frac- 
tion /icM is sensitive to physical treatment, varying 
by 40% over t he range of treatmen ts studied. The 
ART models of iKravtsov et all ()2005D exhibit low val- 
ues /icM = 0.58 ± 0.08, but most (8 out of 9) of their 
halos are group-scale systems with kT < 3 keV. 

Given the degree of uncertainty in the current genera- 
tion of simulations, we quote only a weak constraint on 
the halo ICM fraction 



(/icmTT) 

S 



0.80 ±0.15. 



(14) 



Note that the la range nearly allows for the hot phase 
baryon fraction to identically match the global value. In 
the outer parts of clusters, the scale dependence of the 
hot gas fraction is weak, and we assume the above to 
hold over the range A ^-^200 — 500. 
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Fig. 11. — Ratio of galaxy to dark matter velocity dis- 
persion within r2oo for galaxies with absolute magnitudes 
Mb < —17 in the catalogues produced for the Millen- 
nium Simulation by De Lucia & Blaizot (2007). Results 
are shown for all halos more massive than 10^^h~^MQ 
at z = 1 (upper panel) and z ~ (lower panel). 

The galaxy velocity bias determined by direct simu- 
lation is a subject of active investigation. By identify- 
ing galaxie s with sub-halos in a high resolution A'^-body 
simulation, IColfn et al.l (|200Cl( ) found by = 1.2 — 1.3, a re- 
sult that conflicted wi th early gas dy namic simulations 
that indicated by < 1. iDiemand et al.. (2004b ) point out 
that subhalo populations in clusters are very close to 
stationary solutions of the Jeans equation. The mas- 
sive sub-halo population with significant velocity bias is 
therefore more extended than the dark matter, and is 
in conflict with the observed number density profiles of 
cl uster galax i es. 

iGao et all (|2004[ ) show that galaxy populations whose 
formation and evolution are followed within high- 
resolution dark matter simulations using the semiana- 
lytic approach pioneered by [Springcl ct al. (200_L) can 
match the observed galaxy density profiles of clusters 
well and, as a result, have smaller kinematic bias. Fig- 
ure 11 gives the ratio of galaxy to dark matter veloc- 
ity dispersion for all halos in the Millennium Simulation 
more massive than IO^/i-^Mq. Results are shown at 
redshifts z = and 1, and for all galaxies with ab- 
solute B-band absolute magnitudes Mb < —17. The 
galaxy catalogues used here are based on the galaxy 
formation model of De Lucia & Blaizot (2007) which 
gives good fits to the observed field and cluster lumi- 
nosity functions and to the masses of observed central 
cluster galaxies, as well as fitting the observed galaxy 
number density profiles of clusters. The mean velocity 
bias changes from being slightly less than unity at high 
redhift, (6„) = 0.959 ± 0.002, to shghtly above one at the 
pres ent epoch, (&^) = 1.037 ± 0.00 1 . 

iFaltenbacher fc Diemandl (|2006[ ) show that an alterna- 
tive approach of mapping galaxy luminosity to sub-halo 



size at the time o f accretion (|Nagai fc Kravtsovl l2005l : 
iConrov et al.ll2006f ) yields a similarly we a k velo city bias, 
by = 1.02, within rvar- iBiviano et al.l lj2006[ ) analyze 
vi rial mass estimates u sing the gas dynamic simulations 
of lBorgani et al.l (|2004f ). and find a small, negative effect, 
6„ = 0.95. 

In short, the latest simulations that use either phys- 
ically motivated sub-halo models or a direct approach 
to baryon cooling and star formation indicate that the 
velocity bias is likely to be small. We summarize the 
current situation as 



{bv)s 



1.00 ±0.05. 



(15) 



4.3. 



Cluster Space Densities, Gas Fractions and the 
Power Spectrum Normalization 

The space density of clusters hotter than 6 keV and ly- 
ing within 2: ==0.09 is observ e d to b e 10"'5-^=^°-^/i^ Mpc"^ 
(|Ikebe et al.l 120021: iHenrvl 12004^ . This observation 
sparked a flurry of activity focused on determining the 
power spectrum normalization (Tg (see, e.g., Table 4 of 
Henry 2004). Over the years, published values tended to 
vary by many times their statistical uncertainties. The 
systematic variation was driven partly by different sta- 
tistical treatments and partly by differences in the a s- 
sumed mass-temperature relation ([Pierpaoli et al.l2003l . 
Mass normalization uncertainty introduces error in the 
ower spectrum normalization, Sjlnas) 0.4(S(lnM) 
Evrard et ai][2?M H^rer fc Whitel[200l . 

We now revisit this issue from the inverse perspective. 
We aim to derive the mass-temperature normalization, 
or the equivalent b^, required to match the observations. 

The halo space density is now well characterized as 
a function of the spatially filtered , linear matter power 
spectrum (ISheth fc TormenI Il999t iJenkins et all l2001l: 



iReed et al.ll2007HWarren et al.ri2006[ ). For massive clus- 
ters, the normalization as and matter density parameter 
Qm control the space density via the combination 



^8 = fT8(r!„/0.3) 



0.35 



(16) 



This feature is demonstrated in Figure [121 Instead of 
mass as the independent variable, we plot the cumulative 
space density as a function of potential well depth , 
determined by convolving the M200 mass function with 
the dark matter virial relation determined above. 

The three line styles show matter densities ftm — O-^i, 
0.30 and 0.36, using mass fu nction parame t ers ex pressed 
as linear functions of Qm (jEvrard et al.l l2002[ ). The 
high normalization, 5*8 = 0.90, is approp riate for a con- 
cordance cosmology (jBahcall et al.|[l999() while the low 
value, ^8 = =0.69, is based on the recent WMAP3-hSDSS 
analysis of lTegmark et al.l ()2006f ). who find f}„i = 0.239 ± 
0.018 and (78 = 0.756 ± 0.035. 

The fixed 5*8 condition sets the velocity dispersion, or 
halo mass, at which the cumulative halo space density 
matches the observed space density of 6 keV clusters. 
Assuming /3spoc = l, the 6 keV threshold is equivalent to 
velocity dispersion (Tgai = 988 km s^^. The filled squares 
in Figure [12] are estimates of the cluster space density as 
a functi on of galaxy velo city dispersion from the CIRS 
sample (jRines et al.ll2006l ). 

The CIRS combines local X-ray-selectcd cluster loca- 
tions with the SDSS spectroscopic galaxy sample, which 
is 93% complete within the SDSS sky area for galaxies 
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Fig. 12. — Cumulative space density of halos as a function of dark matter velocity dispersion within r200i calculated at 2 = 0.06 for 
two values of Ss = crg{Qrn/0.3y^''^^ listed. The predictions res ult from convolving the virial scaling, equation JgJ and Table |3] with the 
Jenkins mass function calibrated for M200 l|Evrard et al.ll2002l ) using model parameters appropriate for r2m = 0.24 (dashed), 0.30 (solid), 
and 0.3 6 (dotted). The shaded region shows the observed, local space density of clusters with X-ray temperature greater than 6 keV 
l|Ikebe e t al. 2002; Henry 20Q1). Solid squares with error bars show the space density of CIRS clusters as a function of galaxy velocity 
dispersion fjRines et al.ll26C)M ). with horizontal error bars showing the effect of a ±5% velocity bias. 
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with r-band magnitudes brighter than r— YI .77 . Velocity 
dispersions are measured within values of r2oo estimated 
directly from the sample virial masses. The space den- 
sity is determined from 72 clusters in a volume-limited 
sample within z = 0.10. Note that the number density 
above 1000 km s~^ agrees well with the 6 keV X-ray 
abundance, implying mean /3spec = l- 

Although there are likely to be additional biases due 
to projection and the specific r-band galaxy selection, we 
do not attempt to model these separately here, and the 
inferred measures of 6u for CIRS must be interpreted 
as incorporating these effects. Recent work on mock 
SPSS cluster sampl es defined by a maxbcg algorithm 
(jKoester et al.l [20071 ) suggests that the line-of-sight dis- 
persion may be 10% low compared to ggai o f the best- 
match contributing halo (jBecker et al.l[2007f). implying 
an eff ective value of 6„ = 0.9. However. iBiviano et al.l 
()2006f ) see a smaller effect using hydrodynamic simula- 
tions of galaxy formation. The implications for CIRS are 
unclear without targeted modeling that tests for sensitiv- 
ity to details such as cluster selection method, choice of 
radial scale, redshift and richness. 

Figure [T^ shows that the CIRS data follow the high 
5*8 predictions if the overall bias in galaxy velocity dis- 
persion is small, 6„ = 1. The low normalization mod- 
els require a substantial velocity bias hy = 1.25. The 
recent direct simulations discussed above do not admit 
such a large effect. The horizontal error bars in Fig- 
ure [12] reflect the estimate of the theoretical uncertainty 
quoted in equation (|15|) . Galaxy velocities need to be 
substantially enhanced to allow the low S'g normaliza- 
tion. An effect of this magnitude a ppears unlikely from 
line-of-sight projec tion effects alone (jBiviano et al.|[2006l : 
iBecker et al't[2007l) . 

For the X-ray observations, the gas thermal energy 
parameter &y serves as the tunable parameter to map 
onto the obse rved temperature fun ction. For the concor- 
dance model, lEvrard et all (|2002[ ) find = 1.08 ± 0.07 
for fT8 = 0.9, with a oc ag,~^/^ scaling driven by the 
behavior of the tail of the mass function. We generalize 
this result to non-concordance models by replacing cts 
with ^g. 

Before comparing this result with the range allowed 
by gas dynamic simulations, we need to apply a correc- 
tion to account for the radial scale of X-ray temperature 
measurements. We have calibrated the virial relation at 
A = 200, but the observed values of Tx are derived from 
photons that sample regions interior to this. The exact 
regions sampled depend on the observing co nditions, but 
high-quality data typically extend to rsoo (jPratt et al.l 
[20061 . 

We thus apply a correction factor to 6|n that aligns the 
dark matter velocity dispersion to rsoo- We use a factor 
f^DM(< ''20o)/crDM(< ^500) = 0.84 ± 0.06, where the cen- 
tral value assumes an NFW profile with concentration^ 
C200 = 4.8 (jPolag et al.ll200^ . and the quoted error re- 
flects varying the concentration between values of 3 and 
15. 

The scale-aligned value of the speciflc thermal energy 

^ Note that the Dolag et al. (2004) concentration values are 
quoted against thresholds defined with respect to the background, 
not critical, density. 



ratio that matches the local temperature function is 

&Tlr5oo = (0.77 ±0.07) S-"\ (17) 

This temperature function constraint equivalently sets 
the zero-point of the mass-temperature relation 

(M5oo(6 keV)) = (0.54 ± 0.05) Sl'^ x 10^^ Mq. 

This result, when combined with X-ray observations 
of the scaling of hot gas mass with temp erature, de- 
termin es the ICM mass fraction at 6 keV. iMohr et al.l 
(|1999( l measure ICM masses at rsoo for 45 nearby clus- 
ters in an X-ray flux-limited sample. Before forming a 
mass fracton, we need to align the A/500 normalization 
of that work to the value used here. The total mass 
normalization assumed by Mohr et al. is Af5oo(6 keV) = 
0. 46x10^^/1^^ M0. Aligning with the mass scale derived 
from the temperature function, equation (|18p . requires 
an increase in mass of 15%, or 5% in linear scale. We 
use the fact that observed surface brightness profiles in- 
dicate a nearly isothermal I CM density profile near rsoo 

eumann fc Arnaudl[200l . so the gas mass will scale 
linearly with radius. We therefore multiply the published 
result by 1.055(T8^/®. Note that recent Chandra analysis 
of eleven relaxed clusters indicates that the gas density 
profile may be somewhat steeper than isothermal, but 
this affects our correction at the one-percent level. 

However, we make a second correction for gas asymme- 
try/clumping that essentially nullifies this 5% increase, 
while retaining the 5*8 scaling. Spherically symmetric 
models provide an upper limit on the gas mass within a 
radial shell, as any asymmetries will reduce the amount 
of gas needed to produce the same emission measure. 
Simulations show a small degree of gas clumping in- 
terior to r2oo that gro ws rapidly beyond that radius 
(iMathiesen et all Il999l ). Recent ART simulations by 
iNagail ( 2006( ) displav an overestimate in gas mass of 5% 
at rsoo. 

Adding 5% additional uncertainty in these corrections, 
the ICM mass normalization at 6 keV is 

(A'/icm(6 keV)) = (0.64 ± 0.05) Sl'^ x lO^^'/ifo^/^ Mq. 

(19) 

A potential source of bias that we do not attempt to 
correct for here is that arising from sample selection. 
It is likely that local, X-ray flux-limited samples are 
missing the low-lu minosity tail of the halo population 
([Stanek et al.ll20d6[ ) , and this may correspond to the sub- 
set with low ICM gas fractions. However, proper model- 
ing of this effect requires that we understand the covari- 
ance between Micm; T^, and for the cluster popula- 
tion. We leave this issue for future work, as it is not a 
limiting factor in the conclusions drawn below. 

With equations (fTS]) and (fT9|l . we now form the mean 
ICM mass fraction at 6 keV based on simultaneously 
matching the observed temperature function and ICM 
mass-temperature relation 



(M] 



icm; 



(M, 



500; 



(0.083 ± 0.009) hj^'"^ 



5/3 



(20) 



We can rewrite this result in terms of the halo baryon 
fraction and ICM mass fraction, equation ([Hj), and then 
apply the WMAP3 constrain t on the global bary on frac- 
tion flb/nm = 0.17 ± 0.01 (jSpergel et al.ll2006f) . This 
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Fig. 13. — Dependence of hot gas thermal energy ratio (upper) and mass fraction (lower) on the power spectrum normalization, derived 
from empirical constraints based on the observed space density of clusters above 6 keV and the observed ICM mass— temperature relation. 
Shaded regions show expectations for these parameters from direct numerical simulation. The vertical hatched region is a combined 
WMAP3+SDSS constraint on erg 

(n„i/0.3)i/2 from lTegmark et al.l i200^ . All regions are shown are 90% confidence intervals. 



constrains the parameter combination 

/i?(^7iCMT(6 keV) = (0.49 ± 0.06) ^s"'/^ (21) 

We plot this constraint, along with that derived for 
the specific energy, equation (|17p . in Figure [131 Bold 
lines give the central result while thin lines show 90%- 
confidence uncertainties. The horizontal, shaded regions 
along these curves show the 90%-confidence parameter 
ranges favored by direct simulations, equations (jl3[) and 
equation (|14|) for the X-ray temperature and gas mass 
fraction, respectively. For the spectroscopic to mass- 
weighted temper ature conver s ion, w e assume (Tx/Tm) = 
1.10 ± 0.05 from lNagai eTall (|2006l) . 

The vertical shaded region shows the 90%-confidence 
range of S-a allowed by the W MAP3 and SPSS lum i- 
nous red galaxy sample analysis of lTegmark et al.l l)2006l ) . 
This estimate favors a low mass normalization for cluster 
halos, implying high values for the gas specific thermal 
energy h\ and high values for the hot gas fraction. In 
fact, the latter must be very close to the cosmic value. 

There is significant tension between the low values of 
h\ seen in simulations and the high values required by 
the WMAP3+SDSS analysis. Even after applying the 
spectroscopic to mass-weighted correction, an increase 
of 10%, the gas retains a slightly lower specific thermal 
energy than the dark matter, = 0.96 ± 0.07. Matching 
the temperature function in the best-fit WMAP3+SDSS 
cosmology requires h\ ~ 1.4 ± 0.2, nearly a factor 1.5 
hotter than current models expect. 

Coupled with the galaxy velocity bias discussed above, 
the low normalization forces the conclusion that current 
direct simulations are missing 50% of the specific energy 
in both the hot gas and galaxy components. It is not 



clear how to modify simulation physics so as to bridge 
this gap. Since the two components have very different 
dynamics and thermodynamic behavior, it is likely that 
two or more independent processes must be invoked. 

A higher normalization alleviates much of this ten- 
sion in energetics. Taking 5*8 = 0.80 would require only 
6y = 1.1, a value that lies within 2a of the current the- 
oretical expectation. Slightly higher central tempera- 
tures, perhaps combined with small errors introduced in 
the scale alignment argument above, could increase h\ 
by such a modest amount. The velocity bias parameter 
would need to be positive, but small, 6„ = 1.05. Again, 
this lies within the reach of current direct simulations. 

4.4. Discussion 

The above evidence argues for a higher power spec- 
trum normalization, so that the predicted cluster space 
densities and, to a lesser extent, baryon fractions can 
match observations without severely taxing direct dy- 
namical simulations. Independent lines of evidence also 
point to a higher power spectrum normalization. Joint 
analysis of SDSS, CMB, supernovae, and Ly-ct for est ob- 
servations lead to 5*8 = 0.78 ± 0.03 (|Seliak et all [2006^. 
Weak gravitational shear measurements i n the CTIO 
lensing survey produce a similar constraint ()Jarvis et al.l 
|2006( ). while combined deep and wide sampl es of the 
CFHT Legacy Survey (jSemboloni et al.l l2006( l imply a 
higher value, 5*8 = 0. 89 ±0.06, that reinforces earher work 
(jMasscv et al."2005|). 

LPederscn & Dahl^ (|2006l ) find 5*8 = 0.88 ± 0.09 from 
weak lensing mass determinations of a local ensemble 
of 30 X-ray-selected clusters with t emperatures near 
8 keV. However, subsequent analysis bv lDahlj ()2006[ ) of a 
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slightly extended sample revises this value substantially 
downward, to fi-s = 0.67t5]:[!5. Uncertainty in mean clus- 
ter masses is also apparent in recent studies of the scaling 
of X-r ay luminosity with total mass. While [Stanek et all 
(|2006l ) favor a high Ss ~ 0.8 - 0.9 coupled with a sub- 
stantial level of intrinsic scatter in mass at fixed lumi- 
nosity, <T\nM — 0.25 — 0.4, an alternative perspective of 
5*8 = 0.7 with essentially no scatter is advocated by 
iReiprichl (|2006( l. It is worth keeping in mind that the 
difference in halo mass between 5*8 = 0.7 and 0.9 is sub- 
stantial, a factor or two at Mpc~'^ space density 
(Figure [HI). 

The WMAPS-t-SDSS normahzation imphes that the 
hot ICM mass fraction within rsoo is 10% lower than 
the cosmic baryo n-to-total ratio, a value consistent with 
previous studies () Allen et al.ll200'S ). Higher normaliza- 
tions require this fraction to be lower. For 5*8 = 0.8, the 

value is (0.71 ± 0.09)h^Q^^ . Recent evidence for lower 
gas fractions comes from interferometric observations of 
the Sunyaev-Zel'dovic h effect in a sample of 38 clusters 
(jLaRoQue et al.|[2b06[ ). from analysis of CMB dist o rtions 
induced by nearby clusters (|Afshordi et all 120051 120061) 
and from C handra X-ray observations of relaxed clus- 
ters (jVikhli nin et al1 l2006[ ). Higher normalizations may 
run into difhculties with gas depletion being too large. A 
concordance model, ^8 = 0.9, implies an ICM fraction of 

0.58hjQ^'^Qb/^m within ^500. Such a low value requires a 
higher star formation efficiency, which would be difhcult 
to reconcile with optical observation, or larger baryon 
losses, perhaps driven by AGN feedback. 

Our analysis of the space density as a function of 
ga laxy velocity dispe rsion reaches the same conclusions 
as iRines et al.l ()2006D but our approaches differ. Their 
analysis is based on virial mass estimates that are as- 
sumed to be unbiased. They conclude that a third-year 
WMAP value of Ss requires a substantial velocity bias 
by = 1.30 ± 0.05, a result that agrees well with our Fig- 
ure [12] The fact that this loop closes supports the as- 
sumption that properly applied virial mass est imators are 
not strongly biased (see iBiviano et al.l (|2006[ ) and refer- 
ences therein). 

— 5 /3 

Finally, we note that the common 5g scaling of both 
the gas temperature and the ICM mass fraction is in ap- 
parent contradiction with expectat ions from simple phys - 
ical models for cluster structure (jOstriker et al.l [20051 ). 
Within a given potential, one would anticipate an anti- 
correlation between Tx and /icm, since gas heated above 
the virial temperature would tend to expand and thereby 
reduce the local gas mass fraction. The scalings derived 
here from statistical arguments must ultimately be recon- 
ciled with physical models based on detailed observations 
of clusters at X-ray, optical and sub-mm wavelengths. 

5. CONCLUSIONS 



We present a precise calibration of the virial scaling for 
massive dark matter halos derived from an ensemble of 
large-scale structure simulations. The large majority 
90%) of halos more massive than 10"'^^ Mq comprise 
a regular population that obeys a power-law virial re- 
lation between one-dimensional velocity dispersion ctdm 
and mass, with self-similar slope 0.3361 ±0.0026 and log- 
normal scatter a\na = 0.0426 ± 0.0015. The value of the 
intercept at h{z)M2oo = lO^^ Mq is 1082.9 ±4.0 km s^^ 
The remaining ~ 10% of the population arc smaller mem- 
bers of merging systems (satellite objects) that display an 
elevated mean and dispersion. Our study uses a very sim- 
ple approach to define satellites, and we anticipate fur- 
ther improvements in understanding this minority pop- 
ulation with a more refined halo treatment. The virial 
relation is consistent across six cosmological codes and, 
within the finite range sampled, appears insensitive to 
assumed cosmology or the presence of a trace amount of 
dissipationless baryons. Future work should address the 
robustness of this result to the inclusion of gas that is 
allowed to cool, dissipate and form stars. 

Combining this relation with a precise determination 
of the mass function enables detailed predictions for the 
space density of halos as a function of their internal spe- 
cific energy. Introducing dimensionless ratios to describe 
baryon component specific energies, we apply these space 
density predictions to observations of galaxy velocity dis- 
persion and X-ray temperature. We find that a low nor- 
malization cosmology, as{flrn/0-3)°'^^ = 0.69, suggested 
by WMAP3+SDSS analysis requires an increase, by a 
factor 1.5, in values of Tx and a^^^ expected from the 
current generation of direct structure formation simula- 
tions. A higher normalization, 58 = 0.80, alleviates this 
tension and implies that hot baryon fractions within rsoo 
that arc consistent with estimates from current X-ray 
and Sunyaev-Zel'dovich observations. Even larger val- 
ues of the power spectrum normalization would require 
substantial loss of baryons within massive clusters. 
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